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1 Introduction 

Observational astronomy in the X-ray and 7-ray bands of the electromagnetic spec- 
trum is very crucial to explore high-energy physical phenomena in the Universe. X-ray 
or 7-ray observations from the ground-based instruments are not possible due to the 
atmospheric attenuation. In the last four decades, huge efforts have been made to- 
wards the development of space-borne X-ray and 7-ray telescopes. At the same time, 
these efforts are also limited by the hostile space environment, particularly in hard 
X-rays and 7-rays. High energy charged particles coming from outer space and from 
the Solar wind become trapped in the Earth's surrounding magnetic field creating ra- 
diation belts around the Earth, known as the Van Allen radiation belt. Though the 
satellites carrying the X-ray and 7-ray detectors are usually placed below the inner 
radiation belt (altitude varying from a few 100 km to 10, 000 km), still there are some 
localized energetic charged particle regions which may result in severe damages to the 
instruments if they are activated while passing through those regions such as the Po- 
lar Regions and the South Atlantic Anomaly (SAA) region. Apart from these trapped 
charged particle regions, there are the cosmic diffuse radiation and the cosmic rays, 
mainly protons and alpha particles. These cosmic-ray particles, depending on the ge- 
omagnetic strength at the altitude and position, enter into the Earth's atmosphere 
and interact with the atoms and molecules resulting in various secondary particles. 
These primary and secondary cosmic radiations and charged particles bombard the 
detector and satellite materials and may produce secondary or higher order particles 
in prompt interactions (through bremsstrahlung, pair creation etc.) or/and may ini- 
tiate the detector activation. These radiations and particles will increase the detector 
noise which may mask the original source signal. The high-energy charge particles can 
pass through detector and space craft and may deposit a line of charge in the detector 
volume. This may then be detected in the same way as energy deposit produced by the 
'real' X-ray radiations. It is thus very essential to make an accurate estimate of the on- 
board background noise before designing any space experiment. The background noise 
in an X-ray instrument is mainly due to the interactions of the cosmic-ray protons, 
the albedo protons and neutrons due to Earth's atmosphere, Cosmic Diffused Gamma- 
Ray Background (CDGRB), secondary gamma rays formed in the detector material, 
its frame structure and with the satellite (Dean, Lei & Knight 1991). In addition, the 
long-term activation of the detector materials by these radiations or particles or by the 
particles in the trapped particle region are also responsible for the background noise. 
The background noise varies over a wide range of energy and actually depends on the 
detection capability of the specific detector. These background noise compete with the 
signals due to the interactions of source photons with the detectors. It is therefore 
important to understand the interactions of these background components with the 
detector material and to remove them while extracting the source signal. 

The RT-2 experiment aboard the CORONAS-PHOTON satellite (Kotov et al. 2008, 
Nandi et al. 2009) consists of 4 payloads: three X-ray detectors (RT-2/S, RT-2/G & 
RT-2/CZT) and one processing electronic device (RT-2/E). Detailed description of all 
the payloads and their functionality are given in Debnath et al. (2010), Kotoch et 
al. (2010), Nandi et al. (2010) and Sreekumar et al. (2010). The Phoswich detectors 
(RT-2/S & RT-2/G) are made of Nal (Tl) and Csl (Na) scintillating crystals. Both the 
Phoswich detectors are sensitive to detect high energy X-rays in the energy range of 
15 keV to ~ 1000 keV. The RT-2/CZT detector is a solid-state imaging device, which 
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consists of CZT and CMOS detectors. Both the detectors are sensitive in the energy 
range of 20 keV to 150 keV . 

In the present work, we concentrate on predicting/comparing the background in 
these detectors during the passage of the space craft through the low-background equa- 
torial region (away from the SAA and polar regions) , to estimate the sensitivity of the 
detector. Hence, we make a detailed simulation of the interaction of primary and sec- 
ondary protons, cosmic diffused gamma rays, secondary gamma rays and secondary 
neutrons with the detector volume as well as the whole structure of the satellite car- 
rying the detectors. 

In this paper, we carry out Monte Carlo (MC) simulations using the GEANT-4 
toolkit and highlight the effects of shielding material in calculating the background 
noise due to cosmic-ray photons on the detectors of the RT-2 Experiment. In the 
next section (§2), we describe the typical spectrum of the CDGRB, primary cosmic- 
ray proton and secondary proton, gamma-ray and neutron spectrum in the low Earth 
orbit of interest and of a standard GRB source. We also discuss the physical processes 
that are involved while the incident particles and radiations interact with the detector 
material. In §3 & §4, we present the simulation results of the RT-2/S (RT-2/G) and 
RT-2/CZT payloads. In §5 we compare the predicted result with an observed data set. 
Finally, we conclude in section §6. 



2 The Simulation Attributes 

2.1 Detector characteristics & specifications 

The RT-2 Experiment aboard the CORONAS-PHOTON satellite is a unique experi- 
ment as it consists of different types of detectors, namely a scintillator detector (Nal, 
Csl crystal), a solid-state detector (CZT) and a photo-diode (CMOS) detector. The 
specifications and materials used for all these detectors are given in Table 1. 

Table 1: Detector specifications of RT-2/S (RT-2/G) and RT-2/CZT payloads. 



Payload 


RT-2/S (RT-2/G) 


RT-2/CZT 


Detector type 


Nat + Cst 


CZT & CMOS 


Material composition 


Nat (Tl activated) 


Cdo.9Zno.1Te 




Csf (Na activated) 


CMOS (Amorphous Si photo-diode 
Gd 2 2 S:Tb) 


Thickness (mm) 


3 + 25 


5 & 3 


Size (mm) 


116 dia 


40x40 & 24.5x24.5 


Filter (mm) 


Al (0.5) 


graphite (1.0) (protective cover) 


Effective area (cm 2 ) 


105.6 


48 (3 modules) & 5.7 
6' - 6° 


FOV 


4° x 4° (6° x 6°) 


Readout 


PMT 


pixels 



For the simulation of the effects of various background components in the detec- 
tors we consider an approximate mass distribution of the whole CORONAS-PHOTON 
satellite hosting various detectors on it. Figure 1 depicts the simplified shape of the 
whole satellite including the detectors on it. This shape was used in the simulation. 
The major contribution of its mass is from the satellite shell structure which is made 
up of Aluminum (Al) and the electronics modules inside it, which consist of Aluminum 
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(Al), Silicon (Si) and Copper (Cu) as the major elements. We distribute these materi- 
als throughout the satellite cavity for simplification. Also for the detector components 
other than RT-2/S, RT-2/G and RT-2/CZT we consider simplified structures consist- 
ing of the approximate weights of the major components of concerned detectors. The 
overall height of the satellite construction under our consideration is ~ 324 cm and ra- 
dius is ~ 85 cm. For the detectors of our concern i.e., RT-2/S, RT-2/G and RT-2/CZT, 
we use more detailed geometry described in the following sections. 




Fig. 1 A simplified 3D view of the approximate mass distribution of the satellite containing 
the detectors used for the simulation. 



2.2 Primary Particle Generation 

We simulate five major components which cause background noise in X-ray or 7-ray de- 
tectors. They are: the CDGRB photons, secondary gamma-ray photons due to Earth's 
atmosphere, primary Cosmic-Ray (CR) protons, secondary protons due to interaction 
of CR in Earth's atmosphere and the albedo neutrons from the Earth's atmosphere. 

To simulate the CDGRB photons, we generate the incident photons from the surface 
of a hemisphere (at +ve z-axis) of radius 120.0 cm. The center of the whole geometry 
is at the center of the satellite common mounting plate on which the RT-2 detectors 
and some other detectors are mounted (see Figure 1). The position of the incident 
photon generating hemisphere is such that it covers the whole region above the Earth's 
horizon. The randomness of the incident photons are ensured by placing their origin 
randomly on the hemisphere and the directions of the photon momenta are also chosen 
to be random within the solid angle subtended by the dimension of the satellite radius 
at the vertex of each incident photon. We are interested in the response of the detector 
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in the incident energy range of 10 keV — 100 MeV. The CDGRB spectrum can be 
represented by the equation (Gruber et al. 1999), 



dN _ 
dE ~ * 



7.877S" 1 ' 29 exp"^ 41 - 13 if E < 60.0 keV 

4.32 xl0- 4 (^)- 6 - 5 + 8.4 xl0- 3 (|j)- 2 - 58 (1) 
+4.8 x 10~ 4 (r§) 2 ' 05 if E > 60.0 keV, 

where, E is incident photon energy and |f is in the unit of counts / cm 2 / s / sr / keV . 
We divide the entire energy range into 500 bins equal in logarithmic scale and simulate 
100,000 photons in each bin to retain a good statistics in the simulation result. The 
incident photon spectrum on the detectors due to CDGRB is shown in Figure 2a. 

For the simulation of the secondary albedo gamma-ray photons, we consider pho- 
tons randomly originated from a hemispherical surface of radius 300 cm with the center 
coinciding with the center of the satellite common mounting plate. In this case, we 
consider the position of the generating hemisphere at the opposite that of the primary 
CDGRB (i.e., at -ve z-axis). The direction of the secondary photons has been achieved 
in the same way as that of the primary CDGRB. In this case also we consider the 
energy range of 10 keV — 100 MeV . Atmospheric gamma-ray line emission, such as 
the 511 keV emission from positron annihilation, has not been considered, while sim- 
ulating for secondary albedo gamma-ray photons. The energy spectrum is represented 
by (Ajello et al. 2008; Mizuno et al. 2004), 



dN _ 
dE ~ * 



1.87x10" 



if E < 200.0 keV 



I E 1-5.0 | I E 

) ~r\33.7j 

1.01 x 10~ 4 (t^v)' 1 ' 3 if 200.0 keV <E< 20.0 MeV ( 2 ) 



7.29 x 10~ 4 ( -j-^y) 2 if E > 20.0 MeV 



in units of counts / cm 2 / s / sr / keV . We divide the entire energy range into 500 bins 
equal in logarithmic scale and in each bin, we inject 100, 000 photons. The incident 
photon spectrum on the detectors due to secondary gamma-ray photons is shown in 
Figure 2b. 

We also simulate the detector response due to the CR and secondary protons. We 
consider the input differential spectra for the downward and upward going protons. 
While the upward going protons are mostly secondaries from the Earth's atmosphere, 
the downward going component also contains the primary CR proton above the cutoff. 
We consider these spectra near the equatorial region. 

For the downward going proton component we consider random protons from a 
spherical section of radius 300 cm around the satellite with its center coinciding with 
the center of the common mounting plate and covering the open region above the 
Earth's horizon. The randomization in the direction of the particles is the same as 
described for the CDGRB or secondary photons. We carry out the simulation for the 
proton energy range of 100 MeV — 20 GeV . To produce the energy distribution we 
consider the spectral data given by Alcaraz et al. (2000) for the low geomagnetic 
latitude (0 < &m < 0.2). We fit this data using the functional form described by 
Mizuno et al. (2004) presented by the equation, 



|f = 1.23xl0-«(^)-%xp-(^)- a+1 

b 



+16 9 x in- 7 ( E + Ze A x {E+Mc2 ^-[ Mc ^ 2 .. .. x 1 
+ lt>.y X 1U ^ OeV J X (E+Mc>+Ze<t>y*-(Mc^ X f 



(3) 
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in units of counts / cm 2 / s / sr / keV , where e is the magnitude of the electron charge, Z is 
the atomic number of the particle, a = 0.155, El cu t = 5.1 x 10 5 keV, Solar modulation 
4> = 6.5 x 10 5 kV (a value near Solar activity minimum), Mc 2 is the proton mass and 
E2 cut = 12.25 x 10 6 keV. In this case, we divide the whole incident energy range in 
100 bins equal in logarithmic scale and inject 100, 000 photons in each bin. Figure 2c 
presents the incident energy spectrum for the Primary CR protons. 

The contribution of the upward going proton spectra from the Earth's atmosphere 
is considered as follows. For the positional and directional aspects of the generation of 
these protons we used the same methods as the albedo photons. Here in this simulation 
we consider the energy range of 100 MeV — 6 GeV . The spectral form is given as 
(Alcaraz et al. 2000; Mizuno et al. 2004), 



^ = 1.23xlO- 8 f^r C xp-m 



-o+l 



dE \Gev) CXp -fcJ (4) 

in units of counts / cm 2 / s / sr / keV , where, a — 0.155 and E cu t = 5.1 x 10 5 keV. The 
incident secondary proton spectrum is depicted in Figure 2d. 

For the simulation of the secondary neutrons due to the interaction of the CR in the 
Earth's atmosphere, we consider the generation of the neutrons from a hemispherical 
surface of radius 300 cm in the same manner as the secondary photon simulation. The 
incident particle direction is also randomized in the same way. The energy range of 
the neutron simulation is 10 keV — 1 GeV. The spectral form of the neutron energy 
distribution is given by (Armstrong et al. 1973), 



dN 
dE 



9.98 x 10~ 8 (cfi/) -0 ' 5 if 10 keV < E < 1 MeV 
3.16 x 10~ 9 (sf^)" 1 ' if 1 MeV < E < 100 MeV (5) 
3.16 x 10~ 10 (cfr/)" 2 if 100 MeV < E < 100 GeV 

in units of counts / cm 2 / s / sr / keV . We divide the whole energy range in 500 energy 
bins equal in log scale and inject 100, 000 photons in each bin. Figure 2e shows the 
incident spectrum for the secondary neutron generation. 

The spectrum of energy deposition on the detectors is calculated by normalizing 
the deposited spectrum in the following way. We have N incident energy bins equal in 
logarithmic scale, where Ei is the width of the ith bin. In each bin we are simulating I s 
(here 100, 000) number of particles from a surface of area Aj and over the solid angle Qj 
subtended by the detector to the vertex of the incident photons. Now the total number 
of photons in each incident energy bin is Ii — J E J A ^ J n ^dE. Then we calculate the 

normalization constant Ci — -f 1 . To calculate the normalized deposition spectrum we 
divide the whole deposition energy range in M bins. We have D^j number of photons in 
the jth bin of the deposition spectrum due to I s photons in the ith incident bin. Then 
we calculate the normalized photon counts in the jth bin of the deposition spectrum 
for the all incident photons as 



N 



Dij x Ci 



D , = > -» ~/ (6) 

where, Ad is the area of the crystal (detector) and Ej is the width of the jth energy 
range. 

In the present work, we are only dealing with the sources for the prompt background 
noise. So we simulate here some of the main sources of the prompt background noise 
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Fig. 2 The incident energy spectrum due to the (a) CDGRB photons, (b) secondary gamma- 
ray photons, (c) downward going protons, (d) upward going protons and (e) albedo neutrons. 



sources. Apart from the prompt background, significant noise will be present due to 
the detector material activation as the satellite is in a polar orbit. However, there are 
considerable uncertainties in estimating the contribution due to long term activation 
to the total background. This is because the relevant package, namely, Cosima based 
on Geant4 is still in the developmental stage and cannot be trusted for predicting 
backgrounds due to activation (Zoglauer et al. 2008; Zoglauer, 2009). In the present 
paper, we have deferred the inclusion of the effects of activation on the detector. This 
can be dealt with in a future work. 
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To simulate the incident spectrum of a Gamma- Ray Burst (GRB), we consider the 
Band spectral (Band et al. 1993) form of a bright GRB (GRB 880725), which follows 
the equation, 

dN _ f^(Toow) aex P(-l|) ME<{a-P)E 

in units of counts /cm 2 /s/keV, where, A = 7.056 is a constant, a = —0.32 is the low 
energy slope, /3 = —5.0 is the high energy slope and Eq — 67.7 keV is the break energy 
(Strohmayer et al. 1998). The incident photons are generated from a plane of area 
equal to that of the collimator of the detectors and is placed 60 cm above the detector 
base plate. We consider parallel photons falling straight into the detector through the 
collimator. We also consider the source location (GRB) at an angle of ~ 50° to the 
on-axis of the payload and its effect on resulting photon distribution after interacting 
with detector materials. The incident energy spectrum (Eqn. 7) of the GRB is shown 
in Figure 3. 



Incident GRB Energy Spectrum 
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Fig. 3 The incident photon spectrum due to a typical Gamma- Ray Burst (GRB). 



2.3 Physical Processes for Primaries as well as Secondaries 

The CORONAS-PHOTON satellite is placed at an altitude of 550 km in an orbit 
which is inclined at 82.5°, which results in the passage of the satellite through high 
energetic charge particle regions of the SAA, North Cap (NC) and South Cap (SC) for 
~ 40% of its orbital time. In these regions, satellite operations will be restricted for the 
protection of the detectors from the high radiation dose. Apart from the high energetic 
charge particles of SAA, NC and SC regions, the CDGRB, the high energy protons 
and other albedo particles and photons will hit the satellite and detector materials 
to produce the background noise in the detectors either directly or by creating local 
spallation background. The hard X-ray solar flares along with some other astrophysical 
sources (e.g,. GRBs) would be detected by the RT-2 detectors in the energy range of 
15 — 150 keV, extendable up to ~ 1000 keV . The underlying physical processes by 
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which photons can interact in any medium depend on particle energy and the material 
properties of that medium. 

For the simulation of the detectors under this circumstances, we use the Geant4 
simulation toolkit version 9.1.p03 and the cross section data version G4EMLOW5.1, 
G4ABLA3.0, G4NDL3.12. For the electromagnetic processes in the simulation we con- 
sider the low energy electromagnetic physics list. For the photon interactions we are 
using the low energy photo-electric effect (activated with Auger electron production), 
low energy Compton effect, low energy Rayleigh Scattering and low energy Gamma ray 
conversion. In the above mentioned energy range, the incident photon produces elec- 
trons as the secondary particles. Electrons lose their energy through the low energy ion- 
ization, multiple scatterings and low energy bremsstrahlung processes. For positrons we 
consider the bremsstrahlung, annihilation, ionization and multiple scattering processes. 
We are using a production cut-off value of f \im for photon, electron and positron which 
is the same for all the materials. This implies a threshold energy in the Aluminum as 
990 eV, 1.1 keV and 1.1 keV respectively for photons, electrons and positrons. For the 
hadronic interactions we are considering the "LHEP" physics list defined in the Geant4. 
In the Geant4 toolkit, above mentioned processes are defined as separate modules in 
the "process" category and this module is included in the "Physics list" (Agostinelli et 
al. 2003). Further information are available at http://geant4.web.cern.ch/geant4/ 

2.4 Shielding material 

To shield the detectors from X-rays and gamma-rays from off-axis sources and diffuse X- 
rays and other cosmic particles, materials with high atomic numbers such as Lead and 
Gold etc. are used. These high-Z elements, however, have a deep dip in their absorption 
coefficient just below their K-shell binding energies. Sometimes, materials with lower 
atomic number like Tin and Copper are used to absorb the characteristic fluorescent 
X-rays. In the present case, however, the primary objective of the experiment is to 
make very sensitive hard X-ray measurement of solar flares (typically below 50 keV) 
and in particular use the Phoswich technique for background reduction. Since the 
Csl detector in the Phoswich configuration has a good sensitivity above 100 keV, an 
alternate objective is to use the detector as an open monitor above 100 keV for gamma 
ray bursts and other bright sources. For these dual objectives, we find that Tantalum 
has a very good absorption property with atomic number 73, a K-shell binding energy 
of 67.4 keV and high density of 16. 65 gcm~ 3 . For the shielding of the detector we have 
chosen the Tantalum shield in a way which minimizes mass as well as maximizes the 
dual objectives of solar X-ray flare spectroscopy and off-axis source observations. 

3 Simulation of RT-2/S and RT-2/G payload 

The background simulation of the RT-2/S (RT-2/G) detector is carried out with the 
virtual detector constructed within the Geant4 toolkit environment with the construc- 
tion parameters given in Table 2. These parameter values are the same as in the onboard 
detector geometry, but to avoid the complexity in the detector geometry construction 
in the simulation we simplified some of the detector parts which are not expected to 
change the simulation results significantly. These simplifications consist of omitting the 
ribs to hold the collimator and the detector walls, simplification of the construction 
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of the Photo-Multiplier Tube (PMT), etc. A 3D view of the RT-2/S detector which 
is considered in the simulation is given in Figure 4. This figure is the direct output 
from the simulation and agrees closely with the CAD design of the actual detector 
construction. We simulate the detector response for the background noise due to five 
major background components and a typical GRB spectrum as has been described 
earlier. In the following subsection, we will discuss the response of the detector for 
these particular background radiation under our consideration and in the subsequent 
subsection we will give the results for the GRB spectrum. 




3.1 Summary and interpretation of the background simulation results 

In this section we present the energy response spectra (obtained from the simulation) 
in the scintillators Nal and Csl due to the various background components described in 
Section 2.2. These results enable us to see the contribution of various individual back- 
ground components to the total background. In Figure 5a, we have shown the energy 
deposition spectrum in the Nal and Csl crystals for the in-flight shielding configuration. 
The Nal spectrum shows a low energy cut-off below ~ 15 keV. Low energy photons are 
absorbed by the 0.5 mm Al filter used for protection of the crystals. The Csl spectrum, 
also shows a few counts in the low energy range less than ~ 40 keV. These counts are 
mainly due to the partial energy depositions of the higher energy photons. In higher 
energy range, the counts are high due to the photons coming through the collimator 
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Table 2: Detector construction parameter specifications of RT-2/S (RT-2/G) payload. 



Detector Purts 


Dimensions 


13ottom plate 


£t tJ * U A ^ O . U A KJ * O / / 1 


Detector Housing 


20.2 X 22.4 X 26.7 cm A (inner dimension) 
wall thickness 0.3 cm 


Top Plate 


23.0 x 23.0 x 1.0 cm a 
hole radius 6.5 cm 


Collimator 


13.25 x 13.25 x 5.8 cm s (inner dimension) 
wall thickness 0.3 cm 


FOV 


0.02 cm Tantalum plate thickness 
0.4 cm gap between two Ta plates (4° X 4° FOV for S) 
0.6 cm gap between two Ta plates (6° x 6° FOV for G) 


Nal Crystal 


radius 5.8 cm, thickness 0.3 cm 


Csl Crystal 


radius 5.8 cm, thickness 2.5 cm 
(gap between two crystals 0.01 cm) 


Optical Coupling 


(Silicon Oxide) radius 3.8 cm thickness 1.25 cm 


PMT 


Upper part : radius 4.16 cm, height 3.826 cm 
Middle part (conic section) : height 1.024 cm 
Lower part : radius 2.94 cm, height 12.2 cm 
thickness of the whole PMT 0.36 cm (Aluminum) 


Shielding 


0.02 cm thick Tantalum strip, height 1.8 cm, 
around the collimator wall below the FOV to shield 
Nal+Csl crystal. Shielding weight ~ 35 g 



and not interacting on the Nal crystal and as well as due to those photons entering into 
the detector other than the collimator part. The peak in the Nal and Csl spectrum 
around 60 keV is due to the Ta K a fluorescent emission. The detector response to the 
albedo Gamma-ray photons due to the Earth's atmosphere is given in Figure 5b. Here 
in the figure we can see the Ta K a and K$ fluorescent emission in the Nal spectrum 
along with the photon annihilation peak which is also visible in the Csl spectrum. 
Figure 5c presents the energy deposition spectrum in the Nal and Csl crystal due to 
the downward going protons at low Earth orbit position. We have shown the detector 
response spectrum in Figure 5d, for the upward going protons due to the interaction 
of the cosmic rays in the Earth's atmosphere. Figure 5e depicts the energy response 
spectrum in the Nal and Csl crystals for the secondary neutron background spectrum. 
From Figure 5(a,b,c,d,e) we can conclude that the most important contributor to the 
prompt background noise are the CDGRB photons, while the secondary gamma-ray 
photons also promote a significant portion of the noise. 

The results from the simulation of RT-2/G are almost identical (except the nor- 
malization in the low energy range) as that of the RT-2/S, since both of the payloads 
have identical configuration except the FOV of the collimator and 2 mm Al shielding 
on the top of the collimator of RT-2/G to have higher energy cut-off (below 25 keV). 
So, the above results for the RT-2/S are also applicable to RT-2/G payload regarding 
the response of the detector. 



3.2 Response of a GRB spectrum in RT-2/S 

The RT-2 detectors are primarily designed for spectroscopic measurements of solar 
flares in the 10—100 keV region. For efficient background rejection, a thick Csl detector 
was used. While designing the experiment, it was realized that above 100 keV, the RT-2 
detectors will act as omni-directional gamma-ray burst detectors. Here we also simulate 
a GRB source and give the results below. 
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Energy Spectrum Due to CDGRB 



Energy Spectrum Due to Secondary Gamma-Rays 




Energy (keV) 



Energy Spectrum Due to Secondary Neutrons 




Energy (keV) 



Fig. 5 The energy deposition spectrum in the Nal and Csl crystals due to various background 
components incident spectrum of (a) CDGRB photon, (b) secondary gamma-ray photons, (c) 
downward going protons, (d) upward going protons in atmosphere, (e) albedo neutrons. 



We compare the response of the detector for a known source over the total noise 
level due to various background components. So we carry out a simulation for the 
incident GRB spectrum shown in Figure 3. Moreover we consider two source positions: 
on-axis (0°) and off-axis (50°) location with respect to the detector axis. 

In Figure 6(a,b) we have shown the energy deposition spectrum in the Nal and Csl 
crystals due to the GRB incident spectrum for the on-axis (0°) (in red colour) and 
off-axis (50°) (in blue colour) source position. 
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Nal Detector Energy Spectrum Due to GRB Csl Detector Energy Spectrum Due to GRB 




Energy (keV) Energy (keV) 



Fig. 6 The energy deposition spectrum in the (a) Nal and (b) Csl crystals due the GRB 
incident spectrum (for on-axis (0°) and off-axis (50°) source position). 



It is observed that the GRB spectrum in Nal is highly absorbed (due to the FOV 
walls and shielding at the collimator wall) in the entire energy range for off-axis position 
of the source (Figure 6a) and we can also see the peak around 60 keV is due to the Ta 
K a fluorescent emission. 

In case of Csl crystal, the photon energy below ~ 100 keV is getting absorbed 
both for off and on-axis source position. For on-axis, the photons below ~ 100 keV is 
mostly detected by Nal crystal and we can also notice a emission peaks around 30 keV 
due to the K a and Ka fluorescence of Iodine, whereas for off-axis the photons below 
~ 100 keV are partially blocked by the shielding material around the lower part of the 
collimator. Detection of source photon below ~ 25 keV is less significant in Csl crystal. 
This feature is clearly seen for both the off and on-axis source position (Figure 6b). 

In Table 3, we present the number of photons from the incident GRB which deposit 
energy in the Nal and Csl detectors. We have subdivided the whole energy range of 
10 — 1000 keV into smaller energy bands. We have given both the results for the on-axis 
and the off-axis source positions and the counts due to the total estimated background 
(fourth column). The background counts within parenthesis is from the prompt sources. 
To obtain the total backgrounds we need to include the contribution from the long- 
term activations in the detector. In Section 5 we showed that this prompt background 
only accounts for up to ~ 40% of the total background in the Nal crystal. So in order 
to compare the photon counts due to the GRB with the total background noise, we have 
multiplied the background due to prompt emission by a factor of 2.5 and obtained the 
total background and presented in in this column. This extra contribution, given that 
our satellite is polar, is reasonable since even for equatorial orbits the uncertainly could 
be a factor of 2 (Zoglauer, 2009). 

From the results of Table 3 we can see that for the spectrum of GRBs we are 
interested, Nal has a quite large source to noise ratio (S/N) in the energy range 10 — 
150 keV for on-axis case. However, for the off-axis case, the S/N does not permit to 
detect the source at all. On the other hand, for the on-axis case, Csl has a quite high 
S/N value in the energy range of 50 — 500 keV and for the off-axis case this energy range 
gets broadened in 20 — 500 keV. So this simulation of the GRB ensures the capability 
of such source detection by the RT-2/S instrument. Based on the orbital and detector 
configurations a conservative estimate suggests a detection of ~ 20 GRBs per year. 
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Table 3: Photon counts in Nal and Csl Crystals in different energy ranges due to the incident 
GRB spectrum for both the on-axis (0°) and off-axis (50°) orientation of the source posi- 
tion and due to the total estimated background spectrum with the simulated total prompt 
background in the bracket (see text for details). 



Energy range 
in (keV) 


on-axis (0 U ) 
GRB 


off-axis (50°) 
GRB 


total bkg. (prompt) 


Nal 


10 - 20 


1348.9 


2.9 


17.0(6.8) 


20 - 50 


4764.7 


21.8 


45.5(18.2) 


50 - 100 


3449.5 


54.3 


37.2(14.9) 


100 - 150 


996.4 


22.3 


16.8(6.7) 


150 - 1000 


296.3 


26.9 


70.3(28.1) 




Csl 


10 - 20 


21.1 


32.0 


25.8(10.3) 


20 - 50 


167.0 


210.3 


118.1(47.2) 


50 - 100 


512.1 


953.3 


474.6(189.8) 


100 - 150 


692.5 


552.1 


297.9(119.2) 


150 - 500 


848.5 


516.6 


533.6(213.4) 


500 - 1000 


4.8 


6.2 


110.3(44.1) 



4 Simulation of RT-2/CZT payload 

The parameters used to construct the virtual detector of RT-2/CZT payload within 
the Geant4 toolkit environment is given in Table 4. These parameter values are same 
as that used in onboard detector design and we applied some simplifications as we 
have already mentioned in section 3. A 3D view of the RT-2/CZT payload along with 
different components (e.g., detectors, CAM (Coded Aperture Mask), FZP (Fresnel 
Zone Plate), collimator etc.) as considered for the simulation is shown in Figure 7. 
This figure is the direct output from the simulation. 




Fig. 7 A schematic 3D view of RT-2/CZT payload with different components (eg. detectors, 
CAM (Coded Aperture Mask), FZP (Fresnel Zone Plate), collimator etc.). 
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Table 4: Payload construction parameter specifications of RT-2/CZT. 



Detector Parts 


Dimensions 


Bottom plate 


25.0 x 23.0 x 0.8 cm* 


Detector Housing 


20.2 X 22.4 X 26.7 cm* (inner dimension) 
wall thickness 0.3 cm 


Top Plate 


23.0 x 23.0 x 1.0 cm* 


External Collimator 


10.7 x 10.7 x 17.8 cm* (inner dimension) 
wall thickness 0.4 cm 


Internal Collimator 


10.7 x 10.7 x 13.0 cm* (inner dimension) 
wall thickness 0.4 cm. 


FOV 


0.1 cm Aluminum shielded by 0.007 cm Tantalum 
two perpendicular walls to divide the collimator in four quadrants. 


CZT Crystal 


4.0 x 4.0 x 0.5 crn* (3 crystals) 
gap between modules 0.25 cm 


CMOS 


2.4 x 2.4 x 0.3 cm* 


CAM 


0.1 cm thick Tantalum CAM on top of two collimator quadrants. 
1 mm Al filter used for one CAM and other one is open to sky. 


Zone Plate 


0.1 cm thick Ta zone plates at the face and end 
of the collimator on one CZT and CMOS crystals 


Shielding 


0.007 cm Tantalum around the collimator wall. 
0.02 cm Tantalum below the top plate. 
0.02 cm Tantalum around the detector housing side plates. 
Shielding weight ~ 951 g. 



4.1 Summary and interpretation of the background simulation results 

In the same way and for the same purpose as we have discussed for Nal and Csl crystals 
in RT-2/S instrument, we have calculated the energy response spectrum for the three 
CZT and one CMOS detectors in RT-2/CZT. In Figure 8a, we have depicted the energy 
deposition spectrum in one of the CZT and CMOS due to the CDGRB spectrum. For 
the CZT spectrum, we can see a larger amount of fluorescent photons (emission peak 
at ~ 60 keV) due to Ta shielding of thickness 0.2 mm on top plate and around the 
side plates of the payload. The simulated results for the other two CZT modules are 
mostly identical and hence not discussed here. The energy response spectrum in CZT 
and CMOS due to the secondary gamma-ray spectrum is shown in Figure 8b. These 
spectra also show more or less the same features as the primary gamma-ray (CDGRB) 
spectrum. Figure 8c presents the energy deposition spectrum in the two crystals for 
the downward going protons, while Figure 8d shows the same for the upward going 
protons. In these two cases we can see a relatively higher energy deposition in the 
CMOS than in the CZT. The energy response of CZT and CMOS for the secondary 
neutron spectrum is given in Figure 8e. Figure 8(a-e) also shows the dominance of the 
CDGRB and secondary photon spectrum in the noise contribution. 

4.2 Response of a GRB spectrum in RT-2/CZT 

We now simulate the detector for a typical GRB spectrum (Eqn. 7) as the incident 
spectrum shown in Figure 3 for two source location as mentioned in section 3.2. 

In Figure 9(a-d), we plot the energy deposition spectrum separately for the three 
CZT modules and CMOS detector both for the on-axis (0°) and off-axis (50°) GRB 
position. Figure 9a shows the energy deposition spectrum in the CZT (CZT1) contain- 
ing a CAM along with an 1.0 mm thick Al filter in front of the collimator facing the 
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Energy Spectrum Due to CDGRB Energy Spectrum Due to Secondary Gamma-Rays 




Energy (keV) Energy (keV) 



Energy Spectrum Due to Downward Protons Energy Spectrum Due to Upward Protons 




Energy (keV) Energy (keV) 



Energy Spectrum Due to Secondary Neutrons 




Energy (keV) 



Fig. 8 The energy deposition spectrum in the CZT and CMOS due to (a) CDGRB photons, 

(b) secondary gamma-ray photons, (c) downward going protons, (d) upward going protons and 

(c) albedo neutrons. 



sky. The spectrum in the CZT (CZT2) containing a CAM at the face of the collimator 
is given in Figure 9b. Figure 9c presents the same in the CZT (CZT3) containing a 
pair of aligned zone plates (FZPs) placed at 32 cm apart in third quadrant of the 
collimator. In Figure 9d we have shown the spectrum of CMOS detector which is also 
having two aligned zone plates (FZPs) placed in the fourth quadrant of collimator. 
It is evident from Figure 9(a-c) that all the CZT spectra are roughly the same irre- 
spective of their different configurations (Nandi et al. 2010) of the collimator for the 
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CZT1 Detector Energy Spectrum Due to GRB CZT2 Detector Energy Spectrum Due to GRB 




Energy (keV) Energy (keV) 



CZT3 Detector Energy Spectrum Due to GRB CMOS Detector Energy Spectrum Due to GRB 




Energy (keV) Energy (keV) 



Fig. 9 The energy deposition spectrum in the CZT and CMOS crystals due the GRB incident 
spectrum for on-axis (0°) and off-axis (50°) source positions, (a) Spectrum in CZT with a CAM 
along with Al sheet, (b) that in CZT with a CAM, (c) that in CZT with a pair of zone plates 
and (d) that in CMOS with a pair of zone plates. 



off-axis (50°) source position though they differ for the on-axis source position (low 
energy photons are absorbed in CZT1 due to 1 mm Al shielding). This is because for 
the off-axis source position most of the photons which are depositing their energy into 
CZT or CMOS are coming through the detector side plates having uniform shielding 
material. Those photons interacting on the mask pattern are less likely to go to the 
CZT and CMOS to deposit their energy. All the CZT spectra show a prominent Ta 
fluorescent peak (around 60 keV), caused by the shielding material. 

In Table 5, we give the number of photon counts for three CZT modules and 
for CMOS detector due to the GRB spectrum for the on-axis (0°) and off-axis (50°) 
source positions and due to total estimated background spectrum along with the prompt 
background noise as discussed in section 3. We have subdivided the total energy range 
of 10 — 1000 keV into several energy bands. These results could be useful while analyzing 
the data from real observation (GRB) with RT-2 instruments. 

From the results in Table 5 it appears that with CZT detector, we can detect 
GRBs quite confidently in the energy range of 10 — 200keV for the on-axis case and 
50 — 200 keV for the off-axis case. 
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Table 5: Photon counts in the three CZTs and CMOS detector in different energy ranges due 
to the GRB spectrum for both on-axis (0°) and off-axis (50°) source position and due to 
total estimated background spectrum with the simulated total prompt background inside the 
bracket (see text for detail). 



Energy range 


CZT1 


CZT2 


CZT3 


CMOS 


(in keV) 


on-axis (0") CiRB 


10 - 20 


182.1 


936.3 


414.3 


242.9 


20 - 50 


1049.0 


1415.0 


617.5 




50 - 100 


939.8 


1009.8 


438.3 


184.4 


100 - 150 


408.5 


422.6 


186.7 




150 - 200 


141.3 


145.5 


64.4 


0.6 


200 - 1000 


64.7 


66.6 


36.9 






off-axis (50 u ) GRB 


10 - 20 


7.0 


7.5 


7.9 


7.2 


20 - 50 


52.8 


51.0 


55.7 




50 - 100 


210.2 


230.8 


238.4 


20.8 


100 - 150 


173.5 


203.5 


207.8 




150 - 200 


98.9 


130.5 


137.6 


0.4 


200 - 1000 


51.7 


68.9 


68.0 






total bkg. (prompt) 


10 - 20 


2.8(1.1) 


2.5(1.0) 


4.3(1.7) 


2.5(1.0) 


20 - 50 


16.2(6.5) 


17.7(7.1) 


16.3(6.5) 




50 - 100 


59.3(23.7) 


56.5(22.6) 


59.3(23.7) 


11.1(4.5) 


100 - 150 


33.2(14.1) 


34.5(13.8) 


35.9(14.4) 




150 - 200 


24.5(9.8) 


21.2(8.5) 


23.2(9.3) 


1.9(0.8) 


200 - 1000 


34.6(13.9) 


35.8(14.3) 


39.1(15.6) 





5 Comparison with Observed Data 

We compare the simulation results for the various background components with the 
real data measured by the detector in its in-flight operation. We consider one set of 
onboard Nal spectral data of RT-2/S near the equatorial region and far from the 
trapped charged particle regions. This data refers to the satellite position when the 
instruments came out of the polar region, the count rates were found to be dropping 
slowly. We noticed that the low energy count rates stabilized faster than the high energy 
count rates. We took the background region near the equatorial region when: (a) the 
low energy (< 100 keV) count rates were steady (better than a percent stability in 10 
minutes) and (b) when the high energy count rates (> 330 keV) were steady by about 
5% in 10 minutes. 

The primary energy range of the instrument for spectroscopic measurements is ~ 
15 — 100 keV in the Nal detector. Hence we have concentrated on predicting background 
around this region. 

The comparison of the data in the energy range of 20 — 110 keV with the simulated 
background components of primary and secondary gamma rays, downward and upward 
going protons (consisting primary CR and secondary protons) and secondary neutrons 
are shown in Figure 10. 

From Figure 10 we observe that the total simulated background noise due to the 
main sources of the prompt background components are not sufficient to explain the 
total measured background flux. This prompt background components are able to 
address only about 40% of the measured background. The rest of the background noise 
is suspected to be due to the long-term detector activation which is to be considered 
in a future work. Among the simulated components the most significant contributor 
(~ 27%) to the background noise in the NaT crystal are the CDGRB photons, albedo 
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70 80 90100 
Energy (keV) 

Fig. 10 Onboard RT-2/S Nal background spectrum along with various simulated background 
components of CDGRB photons, secondary photons, downward going protons, upward going 
protons and albedo neutrons. 



photons contribute ~ 10%, downward going protons ~ 0.5% and rest from the other 
factors. 

For onboard calibration of the detectors (RT-2/S and RT-2/G), we used radioactive 
source Co-57 (122 keV) that is placed into one of the slats of the collimator. Extensive 
calibration of both the Phoswich detectors (Nal/Csl) on the ground reveals that the 
~ 58 keV emission line (Debnath et al. 2010) is the intrinsic detector background 
feature (not due to the Tantalum (Ta) shielding around the collimator mesh). The 
emission feature of ~ 58 keV, earlier reported by Gruber et al. (1996) was also observed 
in the onboard data of RT-2/S and RT-2/G (Nandi et al. 2009). The extra emission 
line feature at around 20 keV (cutoff energy of RT-2/S is around 15 keV) could be 
due to the detector's electronic noise, which was not observed during ground testing 
and the investigation is still ongoing. These two emission features characteristics will 
be taken up in a separate paper, while considering the long-term activation aspects 
of the background components in the detector. Also the comparison of the simulated 
result with the real data from the other active volumes are under investigation due to 
the calibration concern of the detectors. 

The Csl detector in the phoswich combination is primarily used for background 
rejection and we use it only for very bright sources like GRBs. To investigate the 
background in Csl for faint sources, one also needs to simulate the exact instrument 
characteristics like the pulse shape for partial energy deposition in Nal and Csl. These 
activities would be taken up in a separate paper. We have, however, looked at the Csl 
spectra and ensured that the background spectral shape broadly agrees with the model 
prediction. 
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6 Discussions and Conclusions 

Background simulations of space-borne payloads is one of the challenging tasks to un- 
derstand the space environment as well as the effect of high energy radiation (photons, 
charged particles, neutrons) on the detectors itself. The effects of the major prompt 
background components like cosmic diffused gamma-ray background, secondary gamma- 
ray photons, primary cosmic-ray protons, secondary protons and albedo neutrons on 
the RT-2 payloads (RT-2/S, RT-2/G and RT-2/CZT) are studied in detail, which 
helped the background calibration and source data extraction from the RT-2 Experi- 
ment. The weight of the material (Ta) that is used for shielding purpose for RT-2/S 
(RT-2/G) and RT-2/CZT payload is around 35 g and 951 g respectively. 

The current work of estimating the prompt background noise covers ~ 40% of 
the measured background. The rest part of the background noise is probably due to 
the long-term activation of the detector materials due to the CR or trapped charged 
particles which is to be estimated in subsequent work (Zoglauer et al. 2008; Zoglauer, 
2009) 

As we have already mentioned that this experiment is primarily designed for the 
spectroscopic measurement of Solar flares, but it is also capable of detecting GRBs. 
Since the energy threshold for these GRB detection is about 100 keV (for off-axis source 
positions) (see, Figure 6b), they will be sensitive for GRBs > W~ 5 ergs/cm 2 . The 
probability of on-axis detection is less than 1 in 1000 and hence GRBs are not expected 
to be detected on-axis. As for the uniqueness of the experiment, if we have a reasonable 
spectral response above 100 keV , it is possible to constrain the spectral parameters 
using the data in conjunction with other contemporaneous spectral measurements. The 
exercise of the present paper is to demonstrate that off-axis response can be handled 
in a reasonable way. 

On 30th January, 2009, the CORONAS-PHOTON satellite was launched success- 
fully and all the RT-2 instruments are functioning to our satisfaction. Already several 
gamma ray bursts (Rao et al. 2009; Chakrabarti et al. 2009abc) and solar flares have 
been detected by the instrument. Detailed reports on the on-board data quality and 
backgrounds would be discussed elsewhere. Detailed results on the observed GRBs and 
solar flares are also being submitted for publication elsewhere. 
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